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Abstract 

The classical stimulated Raman scattering system describing reso- 
nant interaction between two electromagnetic waves and a fast relax- 
ing medium wave is studied by introducting a systematic perturbation 
approach in powers of the relaxation time. We separate amplitude and 
phase effects for these complex fields. The analysis of the former shows 
the existence of a stagnation distance after which monotonous energy 
transfer begins from one electromagnetic wave to the other and this 
quantity is calculated. Concerning phase effects we give the condi- 
tions for the formation of a Raman spike from an initial fast and large 
phase jump in one of the waves. The spike evolution and width esti- 
mated from the reduced model agree with the results from numerical 
simulations of the original system. 

PACS: 42.65.Dr Stimulated Raman scattering, 42.55.-f Lasers, 42.65.-k 
Nonlinear optics 



1 



1 Introduction 



The investigation of the transient effects in stimulated Raman scattering 
(SRS) has a long history. The initial interest came from the study of Druhl 
et al in] who examined the interaction of a strong laser pump pulse with a high 
pressure gas. They found that pulses of long duration get completely depleted 
except in some cases where a very rapid and short lived restoration occurs. 
Numerical simulations then showed that this "spike" of pump radiation was 
connected to a rapid vr phase shift or phase flip introduced in the initial 
Stokes wave fTj. 

The theory of SRS is based on a very simple system of equations re- 
sulting from the Pluczek model of interaction of electromagnetic waves with 
molecular vibrations [HI E] This is essentially the wave equation coupled to 
an oscillator equation. Using a slowly varying envelope approximation for 
the fields and assuming the Brillouin resonance conditions on energy and 
momentum, one can obtain the system of three wave interaction between 
two electromagnetic waves and a medium wave (See, for example JI]). This 
system is interesting for its mathematical properties and can be found in 
many fields. First, in the limit of zero damping Chu et al ^U] showed that 
it can be written as zero- curvature conditions for two differential operators, 
so that an inverse scattering transform (1ST) method [12] could be used to 
solve it. However 1ST methods are generally used to solve Cauchy problems 
where one is given all the field data for all space coordinates at an initial 
time to. The theory of SRS leads to a boundary value problem where the 
fields are given for all times at a fixed location xq. As illustrated in [21 H] 
the 1ST method can be developed in this well but one needs to check 

independently that the potential decays sufficiently fast at infinity. This ap- 
proach was recently validated for the model with group velocity dispersion 
by direct comparison with the numerical solution of the problem [J]. 

From another point of view the SRS is an example in optics of the very 
important phenomenon of three wave interaction. A second well known ex- 
ample is the Mandelshctam-Brillouin scattering, where an optical wave inter- 
acts with acoustic phonons. As well as SRS this interaction is used to amplify 
directly laser pulses propagating in an optical fiber ^3], The second har- 
monic generation in non coUinear beams is another example of three wave 
interaction. Overall this phenomenon is present in much of condensed matter 
physics, from the scattering of spin waves on phonons to the nonlinear optics 
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of Bose-Einstein condensates [in]. Therefore the results obtained for SRS 
could be transposed to these other contexts. 

Here we are interested in the limit where the damping of the medium vi- 
brations is large, so that a perturbation theory derived from the 1ST method 
is not applicable. This is precisely when one observes the Raman spike, 
where the pump radiation after being depleted displays a sharp burst of en- 
ergy, which lasts a small fraction of its initial duration. This phenomenon 
first observed in gases can be studied via numerical simulations of the full 
equations [HEj. These show that the Raman spike can be seen only for large 
dissipation, that it appears with some delay and that its amplitude decays for 
increasing propagation distances in the medium. Our aim is to understand 
how this occurs, more precisely 

1. understand quantitatively the energy transfer from one component of 
the field (the pump) to the other (Stokes wave), in particular estimate 
the distance necessary for a given transfer. 

2. explain the time delay of the Raman spike (stagnation effect) 

3. explain the position, the width and the amplitude of the Raman spike 
as a function of the initial phase flip in the Stokes wave 

4. understand the evolution of the fronts of the pump pulse 

We consider a simplified situation where the pump and Stokes pulse are 
not gaussian but almost rectangular. In this case we will see that the effects 
of phase and amplitude can be separated so that the magnitude of the wave 
intensities can be assumed to be independent of time in the spike region, lo- 
cated in the middle of the pulses. The medium variable can be expressed as 
an integral of the product of the amplitudes of the interacting waves. This 
can be expanded using as a small parameter l/F where F is the damping 
of the medium variable. If we limit ourselves to the leading term in this 
series, we obtain the well-known equations describing SRS for continuous 
waves and within this approximation we find the Raman transfer coefficient. 
At the next order we obtain a system of equations where the evolution of 
the pump amplitude is independent from the one of the phase. Though it 
cannot take into account the phase flip effects causing the Raman spike, it 
should be adequate for the transient phenomena occurring on the fronts of 
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the interacting pulses. To get coupling between the phase and amplitudes 
of the pulses we need to consider the third order term in the expansion of 
the medium variable. In this approximation the dispersion of the response is 
taken into account and this causes phase variations to influence the ampli- 
tudes. Thus, we have the simplest approximate theory to describe the Raman 
spike phenomenon. We will show that nevertheless our analytical description 
of the Raman spike formation reproduces the stagnation effect and yields an 
approximate formula for the spike width. The comparison of our analytical 
results with the results of the numerical simulation shows qualitative and in 
some cases quantitative agreement. 

The paper is organized as follows. After introducing the original model 
in section 2, we describe the approximation used and obtain the zero order 
(standard Raman model) in section 3, the front effects corresponding to 
the first order are described in section 4. Dispersion effects corresponding 
to the third order are given in section 5 and analyzed in section 6. We 
present numerical results confirming this approach in section 7 and conclude 
in section 8. 



2 The original model 

The system of equations describing the non-stationary stimulated Raman 
scattering can be written as [12] 

da db ^ 

^ = go, ^ = -Q (1) 
ox ox 

^ + Tq = -gab\ (2) 

where a, b are the amplitudes of the electromagnetic field and q is the ampli- 
tude of the medium. F is the damping coefficient and g is the amplification 
coefficient. This system of equations may be presented in a dimensionless 
form by rescaling a,b and taking as unit of time the dissipation time l/F, 
a = a/g^/^, b = b/g^l'^ and t = t/V. 

One gets, omitting the tildes 

da db , 

^ = go, ^ = -g «, (3) 

ox ox 
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-^ + q = -eah* (4) 

where e = l/F is a small parameter because we consider the case of a large 
dissipation. 

The initial conditions for the system are 

X = 0, a = ao(t), h = ao(t)pexp i(f){t), 
t = 0,q = 

where the last equation indicates that the medium is initially in the ground 
state. The quantity p = e~"' is the initial amplitude difference between a and 
b, in most cases we will assume p << 1 and |a|i;=o to be of order 1 so that 
|a|^=o >> |&|x'=o- The term is a time dependent initial phase. In the 
following we will assume that ao,7 and (p are real. 

Equations and (jH) are such that (t) = \a\'^ + \b\'^ is independent of 
X and therefore D'^{t) = |ap + |6p = aQ(t)(l + p^). There is no symmetry 
in exchanging a and b. In fact the evolution of the field variable q can be 
written 

qxt + qx = eg(|a|^ - 

Even if |ap — |6p|a.=o < there will be no flow of energy from b to a. Note 
also that an initial homogeneous phase in b plays no role in the dynamics. 

We discretized flll2|) in both x and t using Heun's order 2 Runge-Kutta 
method and advance via the following algorithm 

1. given a{x, t) and b{x, t) for a given x compute q{x, t) by integrating ^ 
in t for the initial datum q{x, t = 0) = 0. 

2. advance to a{x + dx, t) and b{x + rfx, t) for all t by integrating (^J) in x. 

3. go to step 1 with x = x + dx. 

The scheme is started at step 1 for x = and the quality of the computation 
is monitored by evaluating the relative error in |ap + In all the runs 
presented it remained smaller than 10~^. 
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3 A systematic approximation: the zero or- 
der 



The equation describing the evolution of the variable q contains the product 
of the amphtudes of the interacting fields. This suggests to use the variables 
cr = ab* and n = a*a — b*b, so that the principal equations take the form 

! = | = 2(,.- + ,V). (5) 

| + (6) 

The formal integration in (jH)) leads to 

/ da d^a d^a \ 

which written as an infinite expansion is an exact solution of (0} that can be 
obtained by integration by parts. 

To see this consider the formal solution of 

q(t) = -e f exp[-(t - t')]a{t')dt' . 

J — oo 

Integrating by parts, this expression can be rewritten as 

q{t) = -e (^a{t) - exp[-(t - t')f-^dt^ , 
or, after a second step, 

da{t) ft . . „,d^a{t') 



q{t) = -eUt)-^ + /_^exp[-(t-t')]- 



dt' 



and so on. 

For t = r we have da{t)/dt = 0{a/r) Thus this expansion becomes a better 
approximation to the solution when the dissipation parameter gets large. 

In the zero order approximation, the medium variable follows instanta- 
neously the field variables 

q = —ea. 
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This is a well known approximation ^3] where the intensities \a\, \b\ can be 
calculated exactly. We get 

m-^i) = D\t)Y^^^^ . (9) 

Here we introduce p = e^"' so that the key quantity in ()8|9j) 

We immediately see that the distance Xc past which a significant transfer of 
energy occurs is 

Xc = ^, (10) 

which we will refer to as the stagnation distance. We will find this quantity 
throughout this work. 

To measure energy transfer between modes a and b, it is convenient to 
introduce the Raman transfer coefficient 

such that < R{x) < 1. Initially a = so that R{x = 0) = 0. As x 
increases a decreases so that R increases towards 1. Assuming ao to be a 
rectangular initial pulse |aoP(t) = H{t — ti) — H{t — where H is the usual 
Heaviside function, we get \aQ\'^dt = t2 — h and 



/+00 
\a\^dt = {t2-ti] 
-oo 



p2g2{l+p2)ex _^ I 



SO that 

2(l+p2),^ _ 

In Fig. [T] we plot |a(x, t)p as a function of time t for x = 3, 4, 5, 6 and 7 
for both the solution of the partial differential equation system ()3I4|) shown 
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Figure 1: Plot of |a(x, t)p as a function of t for different values of x for 
the numerical solution (full line) and the zero order approximation (jH)). The 
parameter 7 = 5. 

in full line and the zero order approximation (jH)) shown in dashed line. The 
initial pulse is given by 

ao(t) = 0.5(tanh((t - U) /wt) - tanh((t - Q M)), (13) 

where U = 20 and te = 80 are respectively the front and back ends of the pulse 
and the front width Wt = 2. The left (resp. right) panel of Fig. Q corresponds 
to 20 < t < 35 (resp. 70 < t < 85). The zero-order approximation predicts 
correctly the field depletion in the center of the pulse but fails for the values 
of t corresponding to the incoming or outgoing fronts. 

Despite this disagreement on the fronts, the estimation (fT^ of the Raman 
transfer coefficient shown in Fig. |21 matches fairly well the one obtained from 
the full numerical solution. Since the pulse is long compared with the typical 
time 1/r (=1 here) the errors made at the front ends are compensated and 
the approximation is good. This would not be the case for a short pulse. 

It seems clear that the transient effects occurring near an incoming or 
outgoing pulse cannot be captured by the zero order approximation and that 



8 



1 



0.8 - 



o 

cc! 



CO 

a3 



0.4 




O 



O 



6 



1 2 



X 



Figure 2: Computation of the Raman transfer coefficient vs. the position 
X for the numerical solution of ()3I4|) . shown as points and the analytical 
expression ()12j) assuming a Heaviside distribution ao(t), shown in full line. 
The parameter 7 = 5. 

higher order terms are needed. In the next section we include these higher 
order terms limiting ourselves to the first order time derivative. 

4 The first order: front eff^ects 

We substitute the corresponding approximation for q in the initial system of 
equation to obtain 



We introduce the real variable p according to the definition p = 2crexp(— 
so that we now only have real variables which satisfy the system of equations 




(14) 




(15) 



dp dp 



(16) 



ox at 




dn 1 dp^ 



(17) 



dx 2 dt 
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dx dt ' ' 

If we remember that a = ah* , one finds that = arga — arg6, and for a; = 0, 
is the function which determines the phase flip process. The system of 
equations given above admits the conserved quantity with respect to x : 

p^ + n^ = D\t)= [al{t) + bl{t)f (19) 

This relation allows us consider only two equations, one for the difference of 
the squares of the interacting waves and another for the phase difference 0. 



5 Analysis of the characteristic equations for 
an incoming front 

We now proceed to study the effect of an incoming front using equations 
()20I2H) . We first simphfy the right hand side of the equation for n by writing 
it n'^ — f{t) where f{t) = D^(D — 2^). Then we write the equations in 
characteristic form 

The equation for n can be integrated using separation of variables and we 
obtain 



and 



n = -^/(t)tanh[^/(t)(x-a;o)], (23) 

t = to+ / ni^d^ = to- log[ = ], (24) 

cosh(y'/(t)xo) 

where the square root can be imaginary so that the itanh becomes a tan. 
The integration constant xq can be computed from the initial condition and 
the relation uq/ D"^ = tanh(7) obtained from the definition of D. We get 

1 

Xo = — atanh[— tanh(7)] . 
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Figure 3: 3D plot of the incoming front of Fig. The left panel shows 
the solution obtained from the characteristics of the approximate model (j22j) 
while the right panel gives the full numerical solution of ()3|4|) . 

The expressions fl23l24j) give n and t in implicit form and therefore need to 
be solved numerically. 

We have integrated numerically the characteristic ordinary differential 
equations (j22|) in x and obtained the solution in the form of characteristics 
t{tQ,x),n{tQ^x). In Fig. |21this estimation is compared to the full numerical 
solution of ()3I4|) for the incoming front of Fig. ^ The results are shown in 
the form of 3D plots |a(a:,t)p, on the left panel is the solution given by the 
characteristics while the full numerical solution is given in the right panel. 
For large x the plots on the left panel are curved because the wave breaks 
i.e. it becomes multivalued. The approximate model (|22j) is therefore not 
sufficient to describe an incoming front. Second derivative terms should be 
included to yield a dynamics similar to the one of Burgers equation. 



6 Dispersion effects: second order approxi- 
mation 

Lets now consider the terms of ((Tj) up to the second order with respect to 
the time derivative and rescale the independent variable x into x' = ex. The 
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substitution of expression 



da d'^a' 



into the equations (jSj) resuhs in the following system of equations: 



da 

dx 



( da d^a' 



(25) 



(26) 



dn 
dx 



— = -2 2|a| 



d\a\ 
~dt 



+ 



d^ 



a\ 



dt^ 



da 



dt 



(27) 



where the primes have been dropped for simplicity. Both these equations 
and the original system of equations (0) and (jHl) have the first integral with 
respect to x rt^ + 4a* a = const{t). 

Once again we introduce the real variable p according to the definition 2a = 
pexp(i0). Using the system of equation (pSjl and ^T7\i one gets the following 
system of equations for the real variables p, n and0 



(28) 

(29) 
(30) 



where again the conservation law (jl9|) holds. 

Since we want to take into account only effects resulting from the phase 
flip process, we will neglect the terms ^ ^ g^^^^-* and Furthermore, we will 
assume that D does not vary with time, i.e. we place ourselves in the central 
part of the rectangular-like initial pulse. Using these approximations the 
system of equations - (jHUj) takes the form: 



dp dp 
dx dt 


d'^p 


= np[l 


dn dp 
dx ^ dt 




fd<P 
[dt 


d(j) dcj) 

— — h n— 

dx dt 


""dt- 




~ dt 




dn dn 2 
dx dt 



dx ~^ dt 



(31) 
(32) 
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These equations can be written in the characteristic form T: 

dt d(f) ^ 

G^i-C doc 



dn 2 
dx 



(33) 



(34) 



From equations (jHHj) it follows that 
characteristics 



= (pito) and the expression to find 

X 

t = to + J niOd^, (35) 



where is the variable of the characteristic. From (j34|) we can compute n(^) 
as 

n{x, to) = -D^ tanh [D^f{h) [x - xo(to)]} , (36) 
where xoito) is a constant of integration and the function /(to) is defined by 



/(«o) 



where the derivative with respect to t is calculated at t = to? and as = 
0(to) = 0(a; = 0,^0)5 its value can be calculated at a; = 0. Thus the function 
/(to) refiects the phase fiip on the initial pulse h of the interacting waves. 
The constant of integration Xo(to) is defined by the initial condition 

n{x = 0, to) =no = D'^ tanh {d"^ f{to)xo{to)} . (37) 

If we use equation (j37|) and the definitions of n and D we obtain: 

no = al{l- exp(-27)) , = al{l + exp{-2j)) , 

then the constant of integration Xo(to) can be rewritten as a;o(to) = 7 [-D^/(to)] ^ 
Substituting this expression into equation leads to 

n{x, to) = -D^ tanh [D'^f{to)x - 7} . (38) 

So equations and (jHHjl allow us to write the expression defining the 
characteristics 



t = tn 



fit 



1 l^gHsh[7-DV(to)a;] 



cosh 7 



(39) 
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Figure 4: 3D plot of the pump amplitude \a\{x, t) obtained by direct integra- 
tion of the original equations (|3l4j) . The initial phase flip is for t = tf = 40 
and the parameters are 7 = 5 and 6 = 0.011. 



From this equation one can extract to as a function of t and x and substitute 
it into to yield the solution of the problem under consideration: 



There is an intrinsic difficulty in this procedure due to the fact that the 
characteristics can cross so that there might not be a unique to ior a given t 
and X. In particular we will see that around the spike region characteristics 
will cross. Nevertheless it is possible to understand what is going on by using 
asympotics in the regions away from where the characteristics cross. We do 
this in the next section. 

7 Analysis of the characteristic equations for 
the Raman spike 

As shown by expressions (jH^ and (PT|) the principal problem is to determine 
the auxiliary function tQ{t,x), from the equation of characteristics (j39|l . To 
do it we will use the numerical solution of this equation. However, it is 



a 



n{x, t) = -D"^ tanh {D^f{to{t, x))x - 7} , 
^{x, t) = 0.5D^ [1 - tanh [D^f{toit, x))x - 7} 



(40) 



(41) 
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possible to obtain some information for small or large x using asymptotic 
expansions of the function tQ{t,x). We write 

t = to + j log[cosh(D2 f + j log[l - tanh{D^ f x) tanh(7)] (42) 

For X small the above expression can be expanded to yield at first order 

t f^to + D"^ X tanh(7) = to + uqx . (43) 

This expression can be immediately deduced from the characteristic equation 
giving dt/dx fl35|) . It is correct even in the spike region as long as f x « 
1. For larger values of x it is necessary to distinguish the region where / << 
leading to a spike from the rest of the domain where / ~ 1. 

Let us consider this case first. Then one can approximate the expression 
(jH} for large x ^2^hy 

t = to + 27tanh(7) - D"^ xtanh(7) . (44) 

Notice how the slope of the t{x) curve, i.e. the direction of the characteristic 
changes as one crosses the stagnation distance Xc = 

The situation in the spike region can be analyzed using the expression for 
n fl40|) . Inside the spike region the argument of the tanh in n, / x — 7 < 
even for large values of x so that there is no depletion. In this region the 
characteristics follow t = to + a; ~ to + ^qx which is very close to for 
large 7. On the sides of the spike < / < 1 so that the argument of the tanh 
can go through but for a larger stagnation distance than Xc- To summarize, 
outside the spike region we expect the characteristic curves to bend smoothly 
from the behavior ()43|) to as x increases while inside the spike region 
they will follow (j43|) . There is also an intermediate region around the spike 
where the characteristic curves will shift to (j44j) but for x > Xc- These sets 
of characteristic curves will cross as we will show in the next section, leading 
to shock formation. 

To understand things quantitatively let us now give a precise form of the 
phase flip function. We chose 

0(x = O,t) = -^ 



1 + tanh 
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where t/ is the location of the phase jump and aj its amphtude. When 
flj = vr we speak of phase flip. With this choice the function /(to) takes the 
form 



The region where the function /(to) has a minimum will be named Raman 
spike (RS) region. If the parameter 6 is chosen to be small, the RS-region is 
very narrow. Outside this region /(to) is equal to one. Hence we can write 



It shows that up to x ~ (7 — 1)/-D^ the value of the intensity of the pump 
wave is a^(a;, t) = and that for x = Xc ~ 7/-D^ this value decreases down 
to O.bD"^. Only inside the RS -region the intensity of the pump wave does not 
vary. This results in the formation of a Raman spike. Fig. 0] shows a 3D 
plot of the pump field |a(x,t)| solution of the original equations ()3I4|) where 
we chose the Heaviside-like initial pump and Stokes distribution (fT^ In this 
case the middle of the pulse is such that D = 1. We took 7 = 5 and clearly 
see in Fig. 2 that the pump field starts getting depleted only past Xc = 5. 
Thus our estimation of the stagnation distance is a good one. 

Let Ars be the width of the RS-pulse given by half of the amplitude, 
i.e. Ars = 2(ti — tm), where the instant ti is defined by the condition 
a^(x, ti) = 0.5D^ and tm is the position of the minimum of the function 
/(to). Using (PT|) and we can obtain the following equation 



Note how the stagnation distance appears again. Furthermore lix <Xc the 
right hand side is negative so that no value of to can be found and there is not 
Raman spike formation. Hence we obtain the stunted growth (stagnation) 
of Raman transformation once again. From the definition of the stagnation 
distance Xc one can see that in dimensional variables Xc is proportional to 
the damping coefficient F and to the initial depletion coefficient 7. 

Let us now proceed to further analyze (jTTj) in order to estimate the width 
of the Raman spike. For that denote 




(45) 



a^(x,t) = O.SL)^ 1 - tanh(D2x - 7) 



(46) 




(47) 
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so that equation (j47p can be written as 



cosh 



6 



(48) 



Since the parameter e is very small because of the smallness of 6, the argument 
of the cosh function in ()48|) is large so that the function can be approximated 
by an exponential. We obtain 



exp 



'to{ti,x) - 1 



f 



6 



so that to (^1; 3;) = t/ + 51og(2/e). If the function to(^; 3;) can be approximated 
by to ~ t — riox., then ti = tf + nQX + 6log{2/e). But in the framework of the 
approximation under consideration tm = tf + n^x, hence Ajis ~ 2(51og(2/e). 
From this expression and the definition of the parameter e one gets 



A 



6 



RS 



log 



Aa} 



X 



(49) 



where it was assumed that 6'^{x — Xc) << Xc, i.e.,. the distance away from 



Xc is small large enough. It is interesting that for x 
RS-pulse approaches the hmiting value 



oo the width of 



A 



RS 



S\ogi2af/5). 



For 6 = 0.5 this limit is equal to 1.27. When the distance x approaches 
the stagnation distance Xc, the width of the RS-pulse increases. Hence the 
Raman spike forms only ii x > Xc, a. feature that we can clearly see in Fig. 

m 



8 Numerical results and discussion 

We present here numerical solutions of the full SRS equations for the spike 
which confirm the analysis done previously. The equations pi2j) have been 
integrated using the method indicated in section 2. We first examine the 
infiuence of the fiip time and amplitude on the formation of the spike. 

Fig. Elshows \a{x, t) p as a function of t for x = 2,4,6 and 8 with a Raman 
spike. The calculation has been done for a fiip amplitude = vr and three 
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Figure 5: Influence of flip time on the formation of the Raman spike. |aP(a;, t) 
for X = 2,4,6 and 8 as a function of t for three different flip times 6 = 0.1 
(full line), 1 (long dash) and 3 (short dash). The amplitude of the flip is vr 
and the other parameters are as in Fig. |3J 

values of flip time 6 = 0.1 (full line) 1 (long dash) and 3 (short dash). As 
expected from the previous section, no manifestation of the phase flip can be 
seen for x < Xc = ~ 4, where Xc is the stagnation distance. Notice how 
the amplitude of the spike remains of order 1 for the fast flip S = 0.1. On 
the contrary for S = 1 the amplitude decays as x increases to 4 and 6. The 
decay is even stronger for 6 = 3 for which the spike is barely noticeable. 

The influence of the phase change on the Raman spike can be seen in Fig. 
El There we computed |a(x,t)p as previously but took 6 = 1 flxed and took 
three values of the phase jump a/ = vr, | and |. The results are similar to 
the ones shown in Fig. No spike can be seen for x < Xc and the spike dies 
off for 6 = ^ whereas it is still weakly present for 5 = f. 

These two flgures show the importance of the ratio af/6. If a f/6 > tt a 
Raman spike is created and subsists for a few stagnation distances. On the 
contrary if af/6 < n the spike is weak and short lived. 

We now compare the position and half width of the Raman spike obtained 
from the direct solution of equations ()lf2|l and the estimates given by the 
characteristic equations. Fig. [3 shows the characteristic curves t{x) for 
different values of to around the spike instant to = tj. As expected from the 
analysis of the previous section, one sees the characteristic curves away from 
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Figure 6: Influence of the amplitude of the phase change on the formation of 
the Raman spike. |ap(a;, t) for x = 2,4,6 and 8 as a function of t for three 
different phase changes of vr (full line), 7r/2 (long dash) and tt/A (short dash) 
occurring during a time 6 = 1. 

the spike region shift smoothly from the behavior ()4H|1 to as a; > Xc- In 
the spike region the characteristics are straight and around the spike one can 
see the transition from (jlHj) to ()44p for a stagnation distance larger than Xc- 
Observe how the characteristics cross for x > Xc indicating a shock. In fact 
there are multiple shocks as shown by the crossings seen for the curves such 
that to < tf in the right top part of the graph. For to > tf all the curves 
seem to accumulate along a line parallel to (03)- This set of curves shows 
the limitation of the our approach which can predict the kinematics of the 
Raman spike as it is created 

n = -D^ tanh[L)2 f{t - x tanh(7)) a; - 7] , (50) 

but breaks down after shocks are formed. In Fig. [7|we plot the maximum 
of the spike observed in the numerical solution of ()ll2j) and one can see that 
this position agrees well with the one given by the accumulation point of the 
characteristics for to > tf. The result is clearly seen in Fig. |S1 where we plot 
the position of the maximum of the spike as a function of x. On the same 
picture we give the t* value such that q{x,t*) = tt/2 and this follows the 
position of the spike. In Fig. El we compare the width of the Raman spike 
with expression ()49|) for a range of values of 6 and see that there is a fairly 
good agreement. 
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Figure 7: Characteristic curves t{x) for different values of to- The parameters 
are the same as for Fig. 0J The curve in dashed hne indicates the position 
of the maximum of the spike in the numerical solution of the full system. 

9 Conclusion 

We have introduced a systematic approximation procedure to simplify the 
stimulated Raman scattering equations where the small parameter is the 
inverse of the damping of the field variable. Since the equations describe 
complex fields, we have studied separately the influence of the amplitude and 
phase of the waves. This provided a better understanding of the phenomenon. 

Concerning the amplitude, the zero order of the approximation yields 
the well-known theory for continuous waves. At the first order we obtain 
a hyperbolic system of equations which exhibits wave breaking. At this 
stage it does not give the accurate front dynamics but it could probably be 
refined to do so. The simple model we introduce allows to understand how 
energy transfer occurs between the different components a and b of the field. 
In particular we show that this transfer is monotonous and occurs after a 
stagnation distance Xc which we specifically compute from the parameters of 
the model. 

For the phase effects we provide a good description of the Raman spike 
phenomenon in the case of waves of constant amplitude through a simple 
system of hyperbolic partial differential equations for the population differ- 
ence n = |ap — |6p and the relative phase of the fields. We show that a spike 
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Figure 8: t position of the maximum of the Raman spike corresponding to 
Fig. 0] as a function of x (full line) . The t position where the phase of q is 
equal to — 7r/2 is in small dashed line. The line t = t/ + x is marked in big 
dash. 

is formed from an initial phase jump of amplitude and duration b in one 
of the fields if the ratio is large enough. From the characteristic curves 
of the system we obtain the time instant and duration of the Raman spike 
as a function of the propagation distance in the medium. These estimations 
are confirmed numerically. 
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